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Based on a sample of 225.3 million J/ip events accumulated with the BESIII detector at the 
BEPCII, the decays of 77' — > tx + tt' -i+r are studied via J ftp —¥ 777'. A clear 77' signal is observed in 
the 7T + 7r _ e + e _ mass spectrum, and the branching fraction is measured to be B(n' — > ir + ir~ e + e~) = 
(2.11 ±0.12 (stat.) ±0.15 (syst.)) x 10~ 3 , which is in good agreement with theoretical predictions 
and the previous measurement, but is determined with much higher precision. No 77' signal is found 
in the tt + tt~ /j + fi~ mass spectrum, and the upper limit is determined to be B(n' — > tt + -k~ fi~) < 
2.9 x 10~ 5 at the 90% confidence level. 

95 PACS numbers: 25.75.Gz, 14.40.Df, 12.38.Mh 



96 I. INTRODUCTION 112 

113 

97 Since the 77' was discovered in 1964 [1, 2], there has 114 

98 been considerable interest in its decay both theoretically 115 

99 and experimentally because of its special role in low en- 116 

100 ergy scale Quantum Chromodynamics (QCD) theory. Its 117 

101 main decay modes, including hadronic and radiative de- 118 

102 cays, have been well measured [3], but the study of 77' 119 

103 anomalous decays is still an open field. 120 

104 Recently, using the radiative decay J/t/j — >• 777' via 121 

105 -0(3686) — s- ir + ir~J/t/j as the source of 77' mesons, 122 

106 CLEO [4] reported the first observation of the conver- 123 

107 sion decay 77' — > 7r + 7r~e + e~, which has been discussed 124 

108 for many years based on the Vector Meson Dominance 125 

109 (VMD) model and Chiral Perturbation Theory [5-7]. 126 
no Theoretically this decay is expected to proceed via ai2? 
m virtual photon intermediate state, 77' — > tt + it~"/* — >i2s 



TT + ir~e + e~, and provides a more stringent test of the 
theories since it involves off-shell photons. In accor- 
dance with theoretical predictions, the two prominent 
features expected for this decay are a peak with a long 
tail just above 2m e in the e + e~ (M e + e -) mass spec- 
trum, and a dominant p° contribution in AI n + 7T - . CLEO 
with limited statistics was unable to explore these dis- 
tributions, although their measured branching fraction, 
6(77' ->■ n+ir-e+e-) = (2.5^ ± 0.5) x 1CT 3 [4], was 
consistent with predicted values around 2 x 10 ~ 3 . In 
addition, the search for 77' — > tt + tt~ fi + fj," , which is pre- 
dicted to be lower by two order of magnitude, was also 
performed. No evident signal was observed, and the up- 
per limit, B(rj ->• 7r + 7r-/i+/j~) < 2.4 x 10~ 4 , at the 90% 
confidence level (C.L.), was determined. 

At BESIII a sample of (225.3±2.8) x 10 6 [8] J/ip events, 
corresponding to 1.2 x 10 6 77' events produced through the 
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129 
130 
131 
132 

133 
134 
135 
136 



radiative decay J/ip — > "/if , was collected in 2009, and of- 179 
fers a unique opportunity to study 77' decays. In addition™ 
to Vj — > tt + tt~1 + 1~ , rj' —> 77r + 7r~ is also studied in orderm 
to determine the ratio of B{rf — > ir + ir~ l + l~) to B{rf —>W2 

B( V '^ir+ir-l+r) m 



77T+7T ). The advantage of measuring g(^<^. 77r + x -) 
is that uncertainties due to the number of J/ip events 
tracking efficiency from tt^ and the radiative photon de- 
tection efficiency cancel. 



II. THE EXPERIMENT AND MONTE CARLO 

190 

SIMULATION 

191 
192 

BEPCII is a double-ring e + e~ collider designed for a 193 
peak luminosity of 10 33 cm _2 s _1 at the center of mass 
energy of 3770 MeV. The cylindrical core of the BE- 
SIII detector consists of a helium-gas-based drift chain- 195 
ber (MDC) for charged track and particle identification 196 
(PID) by dE/dx, a plastic scintillator time-of-flight sys- 197 
tern (TOF), and a 6240-crystal CsI(Tl) Electromagnetic 198 
Calorimeter (EMC) for electron identification and pho- 199 
ton detection. These components are all enclosed in a su- 200 
perconducting solcnoidal magnet providing a 1.0-T mag- 201 
netic field. The solenoid is supported by an octagonal 202 
flux-return yoke with resistive-plate-counter muon detec- 
tor modules (MU) interleaved with steel. The geometri- 
cal acceptance for charged tracks and photons is 93% of 
47r, and the resolutions for charged track momentum and 
photon energy at 1 GeV are 0.5% and 2.5%, respectively. 
More details on the features and capabilities of BESIII 
are provided in Ref. [9]. 

The estimation of backgrounds and the determinations 
of detection efficiencies are performed through Monte 
Carlo (MC) simulations. The BESIII detector is mod- 
eled with the GEANT4 [10, 11]. The production of the 
J/ip resonance is implemented with MC event genera- 
tor KKMC [12, 13], while the decays are performed with 
evtgen [14]. The possible hadronic backgrounds are 
studied using a sample of J/ip inclusive events in which 
the known decays of the J/ip are modeled with branch- 
ing fractions being set to the world average values in 
PDG [3], while the unknown decays are generated with 
the LUNDCHARM model [15]. For if -> tt + tt~1 + 1~ decays, 
a model [16] based on theoretical calculations using the 203 
vector meson dominant model with infinite-width correc- 204 
tions and pseudoscalar meson mixing [7] was developed. 205 



III. ANALYSIS 



208 
209 
210 
211 



The final state in this analysis is 77r + 7r~Z+^~, with2i3 
I being an electron or a muon. The charged tracks imu 
the polar angle range cos 6>| < 0.93 are reconstructed2i5 
from hits in the MDC. Good charged tracks are required2ie 
to pass within ±10 cm of the interaction point in thc2i7 



beam direction and ±1 cm in the plane perpendicular 
to the beam. Photon candidates are reconstructed by 
clustering the EMC crystal energies. The minimum en- 
ergy is 25 MeV for barrel showers (| cos 9\ < 0.8) and 
50 MeV for end-cap showers (0.86 < |cos0| < 0.92). 
To eliminate the showers from charged particles, a pho- 
ton must be separated by at least 15° from any good 
charged track. An EMC timing requirement is used 
to suppress noise and energy deposits unrelated to the 
event. Candidate events are required to contain exactly 
four good charged tracks with zero net charge and at 
least one good photon. To determine the species of the 
final state particles and select the best photon when ad- 
ditional photons are found in an event, the combination 



with the minimum value of y 

1 

xlc 



-1+1- 



is retained. Here 



X'n+n-l+l- = *4C + £j=l XpIdC?) is the SUm ° f the Chi - 

square from the four-constraint (4C) kinematic fit, and 
that from PID, formed by combining TOF and dE/dx in- 
formation of each charged track for each particle hypoth- 
esis (pion, electron, or muon). Events with xlc < ^5 are 
kept as jtt + tt~1 + 1~ candidates. A 4C kinematic fit un- 
der the hypothesis of j2(tt + tt~) is also performed, and 
X 7 2(7r+rr-) > X^+n- l+l- is required to reject possible 
background events from J/ip — > 72(7r + 7r~). 
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FIG. 1: Kinematical distributions for the rf to 7r + 7r - e + e - 
decay: The invariant mass distributions of (a) 7r + 7r - e + e - and 
(b) 7T + 7r~ . Dots with error bars represent the data; the shaded 
area is MC signal shape, the dashed histogram is the rj — > 
7p (vT + 7r _ ) MC line shape, and the solid histogram is the sum 
of MC signal and MC background from rj — > ^p°(ty + tt~). 
Both of these MC simulations are normalized to the yields 
found in Table I. 

A very clear if signal is observed in the 7r + 7r~e + e~ 
invariant mass distribution, shown in Fig. 1(a) after 
the above event selection. MC study shows that the 
dominant background events come from J/ip — > 777', 
rj — > 77r + 7r~ with the rj photon subsequently converted 
into an electron-positron pair; this background is dis- 
played as the dashed histogram in Fig. 1(a). The di-pion 
invariant mass distribution, which is shown in Fig. 1(b), 
shows good agreement between data and MC simula- 
tion. Figure 2 displays the e + e~ mass spectrum after 
requiring |M(7r + 7r- e + e-) - 771(77')! < 0-02 GeV/c 2 ; the 
background from 77r + 7r _ conversions can be easily distin- 
guished. The enhancement close to e + e~ mass threshold 
corresponds to the signal from the 77' — > tt + tt~ e + e~ de- 
cay, and the clear peak around 0.015 GeV/c 2 comes from 
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FIG. 2: The e + e~ invariant mass spectrum of data (dots 
with error bars) after all selection criteria are applied. The 
solid line represents the fit result, the dotted histogram is 
the MC signal shape and the shaded histogram is background 
obtained from 7/ sideband events. 



FIG. 3: The n ir~ fj, ii~ invariant mass distributions of data 
and MC simulation with all selection criteria applied. Dots 
with error bars represent the data, the solid histogram is MC 
signal, and the dashed line indicates inclusive MC. 



255 sponding to the number of events at 90% of the integral 



the background events of rj' — > ^it + it~ where the photomse 
undergoes conversion to an e + e~ pair and the electron257 
(positron) 's momentum is improperly reconstructed as-258 
suming that all the charged tracks are from the inter-259 
action point. The background contributions of J/ijj — >2to 
7r + 7T~7r° and J/ip —> 77r + 7r _ 7r° are estimated from thc26i 



77' sideband region (0.88 GcV/c 2 < 



A/(7T+7r 



22 5 0.90 GcV/c 2 or 1.02 GeV/c 2 < M^+TT'e+e-) 



I < 
1.04 



events, a maximum 



GeV/c 2 ). 

To extract the 77' 
likelihood fit is performed on the observed e + e~ invari 
ant mass distribution with the signal shape described by 
the MC generator specifically developed for this analy- 
sis, the dominant background shape parameterized by a 
smooth function describing the 7 conversion events from 
?/ — > 77r + 7r~, and the contribution (17 events) obtained 
from 77' sideband fixed in the fit to account for the non-7/ 
background. The fit, shown in Fig. 2, yields 429 ± 24 
Tr + n~e + e~ events, and the detection efficiency obtained 
from MC simulation is (16.94 ±0.08)%; both are summa- 
rized in Table I. 

Figure 3 shows the 7r + 7r~ /i + /a - invariant mass spec- 
trum for candidates surviving all selection criteria. The 
contribution from background events, mainly coming 
from J/ip — > ir° ir + ir~~ ir + ir~ and J/ip — > ~/ir + ir~ir + ir~ 
and estimated with the inclusive MC J/ip events, is 
shown as the dashed histogram. Although a few events 
accumulate in the 77' mass region, they are not significant. 

To determine the upper limit on the 77' signal, a series263 
of unbinncd maximum likelihood fits is performed to thc264 
mass spectrum of tt + tt~ fi + fj,~ with an expected 77' signal. 255 
In the fit, the line shape of the 77' signal is determined by 266 
MC simulation, and the background is represented with a267 
second-order Chebychev polynomial. The likelihood dis-268 
tributions of the fit are taken as the probability density269 
function (PDF) directly. The upper limit on the number27o 
of signal events at the 90% C.L. is defined as N v - L , corre-271 



of the PDF. The fit-related uncertainties on 7V U L are 
estimated by using different fit ranges and different or- 
ders of the background polynomial. The maximum one, 
7V U L = 12, and the detection efficiency from MC simu- 
lation, (35.47 ± 0.11)%, are used to evaluate the upper 
limit on the branching fraction. 
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FIG. 4: Scatter plot of M(77r + 7r ) versus M(tt + 7t ) for data. 

The final state is 777r + 7r~ for this mode. The charged 
track and good photon selection are the same as those 
described above, but no PID is applied in the event se- 
lection. A 4C kinematic fit is performed under the hy- 
pothesis of J/ip — > 7r + 7r~77, and xlc < ^ s required. 
For events with more than two photon candidates, the 
combination with the minimum xic ^ s retained. To re- 
ject background events with ir° in the final state, the 
invariant mass of the two photons is required to satisfy 
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FIG. 5: The r yn + ir~ invariant mass spectrum for data after all 
selection criteria are applied. The solid curve is the fit result, 
and the dashed line represents the background polynomial. 



TABLE I: Numbers used in the branching fraction calcula- 
tions: the fitted signal yields, N (or 90% C.L. upper limit); 
the detection efficiency, e. 



rj' decay mode 


e(%) 


N 


TY + Ti~e + e~ 


16.94 ± 0.08 


429 ± 24 


7T + 7T~ jl + jl~ 


35.47 ± 0.11 


< 12 


7 ( o°(7r + 7r _ ) 


45.39 ± 0.07 


158916 ± 425 



304 
305 

272 M(77) > 0.16 GeV/c 2 ; this removes 94% background 306 

273 while the efficiency loss is only 0.73%. The experimental 

274 signature of J/tp 77/ (77' — > 77r + 7r~) is given by the 

275 radiative photon from J/ip decays, that carries a unique 

276 energy of 1.4 GeV. Consequently it is easy to distinguish 

277 this photon from those from 77' decays. In this analysis, 
27s the combination of 77r + 7r~ invariant mass closest to the 

279 rj' mass is chosen to reconstruct the rj' . 

280 Figure 4 shows the scatter plot of M(^n + ii~) versus3i4 

281 A/(7r + 7r _ ) for the candidate events, where the distinct3is 

282 7/ — p° band corresponds to the decay rj' — > 77r + 7r~. A316 

283 very clean 77' peak is observed in the M(77r + 7r _ ) distri-317 

284 bution, as displayed in Fig. 5. The peak is fitted withsis 

285 the MC simulated signal shape convolved with a Gaus-319 

286 sian mass resolution function to account for the difference32o 

287 in mass resolution between data and MC simulations, 321 

288 plus a second-order Chebychev polynomial background322 

289 shape. The fit, shown as the smooth curve in Fig. 5323 

290 gives 158916 ± 425 rf — > jir + ir~ events, and the detec-324 

291 tion efficiency, (45.39 ± 0.07)%, is obtained from the MC325 

292 simulation; these arc tabulated in Table I. In the simu-326 

293 lation of rf — > 7vT + 7r~, since the resonant contribution327 

294 from p° — > 7r + 7r~ is insufficient to describe the data, the32s 

295 non-resonant contribution (known as the "box anomaly" )329 

296 is also included using a decay rate formula [17] deduced33o 

297 from the ones used in Refs. [18-20]. With the parameters33i 

298 tuned with data, the comparison of the simulated dipion332 

299 mass spectrum to data in Fig. 6 shows good agreement. 333 




FIG. 6: The comparison of the simulated ir tt mass spec- 
trum with data. Dots with error bars are data within the 
rj' region ( [0.938, 0.978] GeV/c 2 ), the dashed histogram is 
background obtained from the rj' sideband, and the solid his- 
togram represents the MC simulation. 



IV. SYSTEMATIC ERRORS 

In the measurement of the ratio of the branching frac- 
tions, the possible systematic error sources and the cor- 
responding contributions are discussed in detail below. 

• Form factor uncertainty. In the MC generator 
used to determine the detection efficiency of 7/ — > 
tt + it~1 + 1~ , the VMD factor defined for the hid- 
den gauge model is introduced to account for the 
contribution from the p° meson. The detection ef- 
ficiency dependence is evaluated by replacing the 
factor above with the modified VMD factors de- 
noted in Rcf. [7]. The maximum change of the 
detection efficiencies is assigned as the systematic 
error, which is listed in Table II. 

• MDC tracking efficiency. Since the systematic er- 
rors for the two charged pions cancel by measuring 
the relative branching fraction of 77' — > tt + tt~1 + 1~ 
and 77' — > 77r + 7T _ , only the systematic error caused 
by the MDC tracking from the leptonic pairs need 
be considered. As the momenta of the two charged 
leptons are quite low, it is difficult to select a 
pure sample from data. In this analysis the MDC 
tracking uncertainty of charged pions at low mo- 
mentum is determined and used to estimate that 
of the leptons by reweighting in accordance with 
their momenta. The data sample of J/tJj — > 777', 
rj' — > 77r + 7r~ is used to evaluate the data-MC dif- 
ference of pions at low momentum and finally the 
MDC tracking uncertainty is estimated to be 2.1% 
for electrons and 1.6% for muons, where the domi- 
nant contribution is from the momentum region be- 
low 200 MeV/c. Therefore 4.2% and 3.2% are taken 
as the systematic errors on the tracking efficiency 
for the channels with e + e~ and p + ji~ , respectively, 
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334 in the final states. 

TABLE II: Impact (in %) of the systematic uncertainties on 

• Photon detection efficiency. The photon dctec- the measured branching fractions 

336 tion efficiency is studied with three independent Sources rj -» 7r + 7r~e + e~ rj -» 7r + 7r~/iV~ 

337 decay modes, tp(2S) -> n+n-J/ip (Jjib -> p°7r°), Form factor uncertainty 0.2 0.3 

ib(2S) -> TT+n-J/i; (J/ib -> J+r) and J/^ -> MDC tracking 4.2 3.2 

A [21]. The results indicate that the difference ^° ton detectlon 1.0 1.0 

PID 4 4 

340 between the detection efficiency of data and MC .„ . . ,. „, 

, ,. . .„. r J , , _,. 4C kinematic fit 2.0 2.0 

Simulation is withm l% for each photon. Since Background uncertainty 0.9 

342 the uncertainty from the radiative photons can- / mass w i nc j ow 2 

343 eel by measuring the relative branching fraction of tv 7T + 7r - 0.5 0.5 

344 r)' -> ir + Tr~l + l~ and rj -j> jir + ir~, l% is taken MC statistics 0.6 0.4 

345 to be the systematic error from the photon in rj Total M6 5l3 

346 decaying into 77r + 7T _ . 

• Particle ID. The study of the particle ID efficiency^ 0.5%, contains the error due to the tt° veto cut 
of the pion is performed using the clean control^ (Af (77) > 0.16 GcV/c 2 ) and the fit-related error. 

349 sample of J/ip —> tt + tt tt , and indicates that the 

350 pion particle ID efficiency for data agrees within l%3<m Except for the systematic uncertainties studied above, 

351 of that of the MC simulation in the pion momentum 3 9i a small uncertainty due to the statistical error of the effi- 

352 region. The particle ID efficiency of the electron 39 2 ciencies in vj — > tt + tt~1 + 1~ and rj' — > 7vt + /t~ is also con- 

353 was checked with radiative Bhabha events, and the 3 93 sidercd; all errors are summarized in Table II. The total 

354 difference between data and MC simulation is founds^ systematic error is the sum of them added in quadrature. 

355 to be 1%. In this analysis, 4% is taken as the sys- 

356 tematic error from the particle ID efficiency of the 

357 four charged tracks in rj decaying into ir + ir~l + l~. 395 V. RESULTS 

35s • Kinematic fit. The clean sample J lib — > dm (d> — > _, .. , ,. . . r „, , . . 

„, „_ , _ n . , f , , ' \. The ratio upper limit of B(rj -> ix^ix / + / to 

359 K~<~K , ri — > tt^it 7T U ) selected without a kmc- „/ , + -\ • 1 1 / j -.1 

. ' ' { . B(rj — > 77T 7r ) is calculated with 

360 matic fit is used to estimate the systematic error 

361 associated with the 4C kinematic fit. The differ- »/■„) , _+_-7+/-\ at /, 

* l\ 1\ I Ij ) 77' — >TT~*~ TT ~ l~ I Tj' — yTT~*~ TT~ l~T~ l~ 

362 ence between data and MC is determined to be , , + , = — -. , 

(0.47 ± 1.45)%, with x 2 < 75. In this paper, 1.9% ^ ~* ^ n > ^'-^r+Wty-^+ir- 

is taken to be the systematic error from the kine- where N n+n _ m _ and AL, _ +7r _ are the observed 

matic fit + | h ^ anal y zed ; deca y s , of 3 ,W "VI 397 events (or the 90% C.L. upper limit) of rj -> Tr+Tr-l+i" 

^ ^ ^ ^ 11 • FOT ^ ^ , V 77T+7T d / and £ _ and £ _ &re 

channel, the 4C kinematic fit uncertainty is esti-^ ^ corresponding detection efficiencies. With the mini- 

36s mated to be less than 0.7% using the control sam- . . BiV-s-7r + 7r~e + e~') 

1 tii , rpi ,1 r 1 • r.,400 bers given m Table I, the ratio J,, , r — is dctcr- 

369 pie J/ip — > pTT. lhus, the error from kinematic fit 6 ' B(7j'->77r+7r-) 

370 is, 2.0%, the sum of them added in quadrature. ™ mined to be (7.2 ± 0.4 (stat.) ± 0.5 (syst.)) x 10~ 3 , where 

402 the first error is the statistical error from A r j; /_ >7r + 7r - ; +;- 

371 • Background uncertainty. Studies have shown that™ and iV^_j. 7W + w -. To calculate the upper limit, the sys- 

the mass resolution of 77r+7r- , as simulated by the 404 tematic error is taken into account by a factor of jt . 

373 MC, is underestimated. To evaluate the systematic ™ , - „ in -4 n ""l* 

_ ' . ^05 Iherefore the upper limit, 1.0 x 10 , on the ratio 

374 effect associated with this, the invariant mass of B , , + _ + 

77T+7T- in the MC sample is smeared with a Gaus-™ b W ^+tt-) is S ivcn at thc 90% confidence level. 

376 sian function, where the width of this Gaussian is 

377 floated in the fit. The change of the result, 0.9%, 

378 is assigned to be the systematic error. 407 VI. SUMMARY 

• rj mass window requirement. Another sourcc 408 xlic measurements of rj -> ir+if-l+r, l± = (e ± ,/j ± ) 
of systematic uncertainty is the requirement on 4M are performed using the sample of 225.3 million J/i/> 
the rj mass window selection \M(tt+iT e+e ) - 410 events co n e cted with the BESIII detector. A clear 
m(r]')\ < 0.02 GeV/c . The uncertainty is stud- 411 signal is observed in the invariant mass spectrum of 

383 ied using a looser requirement of 0.90 GeV/c < , _ , _ . d . s(ri '^tt + tt~ e+ e~ 1 . , . , 

_ % i nn n xt 1 1 i 1 • i 412 t 71 e e 1 an d the ratio V, , v — , \ ' is determined 

384 M(7t + tt e + e )< 1.02 GeV/c 2 , and an uncertainty B^'^rr+Tr ) 

of 2.0% is assigned for this item. 413 to be ( 7 - 2 ± 0.4 (stat.) ± 0.5 (syst.)) xJO- 3 . Using 

414 the PDG world average of B(rj — > jtt + tt ) and its un- 

386 • Uncertainty of the number of vj —r jtt + 7t~ cvcnts4is certainty [3], the branching fraction is measured to be 

387 (A^ T? /_ >77r + 7r - ). The uncertainty from this itcm,4i6 B(rj — > 7r + vT _ e + e _ ) = (2. 11±0. 12 (stat. )±0. 15 (syst.)) x 
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417 10~ 3 which is consistent with the theoretical predictions435 

418 and previous measurement, but with the precision im-436 

419 proved significantly. The mass spectra of ir + ir~ and437 

420 e + e~ are also consistent with the theoretical predictions438 

421 that Mj+ T - is dominated by p°, and M e + e - has a peak439 

422 just above 2m e with a long tail. No evidence for r/'tw 

423 decaying into 7r + 7r~/i + /i~ is found, and an upper limit44i 

424 of 1.0 x 1CT 4 on the ratio of g( l^ 7r is ob- 442 

B(ri' — ¥~fir+-n j 

425 tained at the 90% confidence level. The corresponding 443 

426 branching fraction upper limit of rf — > n + n~ /i + pT is 

427 B(if -> n+n-n+fj,-) < 2.9 x 10- 5 . 

446 
447 
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